Quantum Coherence of Discrete Kink Solitons in Ion Traps 
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We propose to realize quantized discrete kinks with cold trapped ions. We show that long-lived 
solitonlike configurations are manifested as deformations of the zigzag structure in the linear Paul 
trap, and are topologically protected in a circular trap with an odd number of ions. We study the 
quantum-mechanical time evolution of a high-frequency, gap separated internal mode of a static 
kink and find long coherence times when the system is cooled to the Doppler limit. The spectral 
properties of the internal modes make them ideally suited for manipulation using current technology. 
This suggests that ion traps can be used to test quantum-mechanical effects with solitons and explore 
ideas for the utilization of the solitonic internal-modes as carriers of quantum information. 



Solitons arc localized configurations of nonlinear sys- 
tems which are nonperturbative and topologically pro- 
tected (H- Quantum- mechanical properties of solitons, 
such as squeezing, have been predicted and measured in 
optical systems Quantum dynamics has been ob- 
served with a single Josephson junction soliton 
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waveguide arrays [3, 01 and Bose-Einstein condensates 
solitons are mean field solutions, localized to a few 
sites of a periodic potential. In chains of coupled parti- 
cles, solitons are discrete spatial configurations, as in the 
Frenkel-Kontorova (FK) model 0,13. 

Discrete solitons of the FK model and its generaliza- 
tions are referred to as kinks. An important property of 
kinks is the existence of localized modes. One mode is the 
kink's translational 'zero- mode', whose frequency gener- 
ally rises above zero. Other localized modes are known as 
'internal modes' 0,0- Physically they describe 'shape- 
change' excitations of the kink and typically they are 
separated by an energy gap from other long-wavelength 
phononic modes. It was suggested to use the internal 
mode as a carrier of quantum information 

Quantum information processing in ion t raps 1 1211 has 



dramatically improved over the last decades 
cently there has been a considerable interest in using 
trapped ions for quantum simulation of various systems 
such as spin-chains 15l4l7l| and Bose-Hubbard models 
field models 19 1, cosmological effects 20 1 and black 
holes 2l|. It was suggested to realize a ID generalization 



of the FK model by adding an external periodic potential 
to an ion trap 22 1 . 

In this Letter we demonstrate that quantum coherence 
in static discrete kinks can be observed with ordinary 
Paul traps without external additions. We explore quasi- 
2D discrete kinks resembling those of the zigzag model 
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FIG. 1: Metastable linear trap configurations with 33 ions, 
(a) An extended kink. The localized internal mode (red bars) 
involves ~ 10 ions, (b) A highly discrete 'odd' kink. A local- 
ized internal oscillation (blue bars) involves ~ 3 ions. 



[23| . In the linear trap we find local metastable defor- 
mations of the zigzag structure 2J|, as depicted in Fig. 
[l] which are long-lived already with a moderate number 
of ions, N > 20. In a circular trap with an odd number 
of ions, similar configurations form the ground state. We 
study the robustness of a high-frequency internal mode of 
the kink against decohercncc in the thermal environment 
of all the other modes. With all nonlinear interactions 
accounted for, we numerically integrate a non-Markovian 
master equation, which leads us to our main result: al- 
ready at the standard Doppler cooling limit coherence 
persists in the internal mode for many oscillations. This 
could allow a first direct measurement of decoherence 
time of solitons. 



Let us consider N ions trapped either in a linear 
trap or in a circular ring trap. Throughout this Let- 
ter we use nondimensional units by employing the nat- 
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ural length scale d = (e^/mi/^) and time scale l/iy, 
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where m is the ions' mass and v is the axial (radial) 
trapping frequency in the linear (circular) trap. In the 
linear case the trapping potential can be expressed as 
V — i {xl + fiyf)- At sufficiently high transverse 
trapping /3, the ions crystalize in a one-dimensional chain 
along the x-axis. When /3 is lowered, the ions undergo a 
second order phase transition [25[ and the lowest energy 
state has a zigzag shape with interesting quantum me- 
chanical properties [2^. However, depending on /3, other 
metastable configurations exist. We have identified two 
additional types of configuration that are not destroyed 
by thermal fluctuations at up to ^ 15 times the Doppler 
temperature. Figure [IJa) shows an extended soliton-like 
deformation (/3 = 40), with a cross-over of the upper ions 
of the zigzag above the lower ones. In Fig [IJb) one ion 
is forced out of the zigzag to form a dense defect at the 
center (/3 = 85). 

In the ring trap, ions are confined at uniform density 
around a circle and the radial confining potential is 



^ 1 



Ri 



N 
^2^ 



(1) 




FIG. 2: Dispersion-relations in the circular trap. At 7 = 0.82 
(red circles and upper inset) the kink is extended, there is no 
gap at the top of the spectrum and the translational mode's 
frequency approaches zero. At 7 = 1.3 (blue diamonds and 
lower inset) the kink is localized on 4 ions and the local- 
ized modes at the top and bottom of the spectrum are gap- 
separated. The mode corresponding to rotational symmetry 
of the entire trap has been omitted. 



where Ri = {xi,yi) and 7 measures the strength of the 
radial trapping, independently of the number of ions. 7 oc 
1/2/3 jg varied by changing the radial trapping frequency. 
The lowest energy configuration with an odd number of 
ions can be one of a few types. At high 7, close to the 
phase-transition from a ID chain (7 « 1.6), the kink is 
localized with two ions facing inside the zigzag shape. 
At lower values of 7 the two kink-core ions lie outside 
of the zigzag, as shown in Fig. [51 lower inset. When 
lowering 7 further, an extended kink is formed (Fig. [2l 
upper inset). The ring topology protects kinks in this 
trap from breaking. 

In the absence of a kink the linearized frequency spec- 
trum of the normal-modes consists of a phonon band ter- 
minated at a cutoff frequency that depends on the ion 
density. The presence of the kink causes a few highly lo- 
calized modes to split away from the rest of the spectrum. 
In a large range of values for the trap parameter, there 
exists a localized mode lying above the band and sepa- 
rated by a gap (Figs. [21 [HJa)). This mode corresponds 
to out-of-phase oscillations of the ions in the kink's core. 
It is this "high-frequency mode" which we suggest for 
coherent manipulations. 

Expanding the Hamiltonian in a perturbative series 
about the classical kink configuration, Ri = R^ + SRi, 
we switch variables to the normal coordinates which di- 
agonalize the harmonic part. Keeping nonlinear terms 
up to the fourth-order we proceed with canonical quan- 



tization of the normal coordinates to get 
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where Hfree is the Hamiltonian of free-phonons, Qj = 
(aj + ftj ), is the creation operator of the normal mode 
J with frequency ujj (a linear combination of all the 5Ri), 
and h ^ 2 X 10"^ is nondimensional. 

We now briefly outline the key steps of a derivation 
of a master equation modelling the coherent quantum- 
mechanical time evolution of the mode of interest. We 



follow the notation in 27| . We divide the modes into the 
'system' - consisting of the high-frequency mode (wi), 
and the 'bath' - all other modes, and split the Hamilto- 
nian H = Hs + Hb + HsB into 3 parts: for the system, 
bath and system-bath interaction, respectively. We next 
apply the Born approximation in the Liouville-von Neu- 
mann equation, assuming the factorization p (0) ® xb 
where O denotes an operator in the interaction picture 
and Xb is a thermal density matrix for the bath, yielding 

P(t) = 

-^tre { [Hsb (t) , Jo dt' [Hsb (f) , p (f) (E> xb]]} (3) 

Expressing the interaction as a sum of terms , each con- 
sisting of a system operator Sa multpilying a bath op- 
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erator Ba, we have Hssit) = fi'^^Sa{t)Ba{t). Assum- 
ing the bath to remain in thermal equihbrium we take 
a free-evolution Hamiltonian for Hb- In Hg, however, 
we keep all terms. In Hsb, we include operators up 
to quadratic order from both the system and the bath. 
Thus, Sa & {9i,6i}, and correspondingly. 
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We define the renormalized bath correlation functions 

C^p (i-t') = ((Sa(i)-(B.>) [Bf,it')~{Bp)))^ (4) 

which do not decay [Fig. [31[b)]. This is due to the dis- 
creteness and cutoff of the bath spectrum, in addition 
to the nonlinearity of the interaction. We therefore pro- 
ceed with a non-Markovian treatment. This results in 
the integro-differential equation 



where Sc^p it) = S„ {t') p (f) Co^p {t - t') dt' . 



(5) 



We solve eq. ([5]) numerically, taking 33 Ca+ ions 
in the configuration of figure [ijb). With v^ jl-K — 
0.88MHz and iyy/27r = 8.1MHz, the mode frequen- 
cies are uji/2tt = 11. 5M Hz for the high-frequency mode, 
W2/27r = 10.6MHz for the next mode, and wiow/27r = 
2.1MHz. The inter-ion separation at the kink center 
is 1.7/iTO. We set the temperature of the bath to the 
Doppler cooling limit TDoppicr = 2tt x lOMHz, and the 



low- frequency mode has 4.3 phonons 28|. In p (0) we as- 
sign a representative superposition state |0) + This 
state can be created after sideband cooling of the inter- 
nal mode and initialization using quantum information 
techniques 13|, [IJI • This must be done slowly compared 
to the inverse of the energy gap separation of this mode 
{2tt X 0.9MHz, or ^ wi/12 in the above example). In 
Fig. |3l[a) we show the fidelity [29^1 of the system's evo- 
lution versus time. The fidelity is calculated [sBl with 
reference to an isolated free phonon, for which there is 
an oscillation of the relative phase between the levels |0) 
and The fidelity remains very high for the simulated 
~ 100 periods of oscillation. This would allow the state 
initialization to be performed with a high fidelity. 

As one test of our results, we compare the master equa- 
tion approach to a direct unitary calculation involving 



only the three modes which are expected to dominate 
the nonlinear process near a resonance as in figure |4l 
The agreement depicted in figure ISjc) indicates that the 
master equation indeed captures the evolution very ac- 
curately, and allows one to simulate the coherent affect 
of the thermal bath modes on the internal-mode. Far 
from the resonances, the Born approximation is valid as 
long as e = A£'i/£'bath ^ 1, where AEi is the energy 
leaking from the system into the bath, and ii^bath can be 
estimated - for the worst case - using the low-frequency 
mode. In the simulation of figure |3lja), e is only a few 
percent, so the Born approximation is justified. 
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FIG. 3: (a) Fidelity of coherent oscillations of the high- 
frequency mode in the linear trap (/3 = 85). (b) i?e(Cii(r)) 
of the same configuration, (c) Fidelity near a resonant con- 
figuration, using the master equation (blue), and a unitary 
evolution (dashed). 



In order to analyze the dependence of the coherence 
on the trap parameter, we take the limit of a large cir- 
cular trap with vanishing curvature, which is also the 
limit of a large linear trap with fixed ion density in the 
center. This is achieved using 'periodic' boundary con- 
ditions such that the longitudinal distance between any 
two ions is evaluated modulo half the chain length. The 
soliton configuration analyzed is of the same type as in 
Fig. [IJb) . Figure IDJb) shows the fidelity as a function of 
7 [which is independent of the number of ions, see Eq. [Tj. 
The fidelity remains high provided that wiow is not too 
low [Fig. mja)], and that there is no strong resonance, 
cji « ujj + wiow , where the largest interaction coefficient 
is with j = 2, a partly-localized phonon. The resonance 
on the left is seen to have a much weaker effect, owing 
to a smaller matrix-element. The loss of fidelity at the 
resonances is mainly due to energy relaxation from the 
high-frequency mode, while away from resonance the off- 
diagonal elements grow considerably, expressing the loss 
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of phase coherence. 
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FIG. 4: Mode frequencies and fidelity analysis, (a) The fre- 
quencies of the lowest and the three highest modes as a func- 
tion of the trap parameter in the 'odd' kink configuration (see 
text for details), (b) The lowest value of fidelity obtained for 
simulations of 50 oscillations of the high-frequency mode, at 
different parameter values. Black diamonds: at inter-ion sep- 
aration of 1.5/im. Blue squares: at twice that distance. The 
numbers indicate the population of thermal phonons in the 
low-frequency mode. 

We continue with a discussion of various aspects of our 
treatment. The geometric coefficients Mijki and succes- 
sive terms omitted from the series expansion in eq. ([2|), 
are not necessarily decreasing compared with the third 
order L.^fe. However, powers of ^ 1/230 dominate 
the convergence at low excitations, so already the con- 
tribution of quartic terms is negligible in general. The 
Hilbert space of the system is truncated at dimensions 
8 — 10, which we have checked to have negligible effect 
on the results because of the low phonon-numbers. 

The dominant contribution to decoherence comes from 
multi-phonon processes involving the localized high and 
low modes. Powers of wiow^^^^ enter the coefficients of 
such interactions which grow stronger as wiow — > 0. Since 
wi depends only weakly on the trap parameter (at a given 
kink type), while wiow and the gap are tunable, the high- 
est fidelity is achieved at the center area of figure SI where 
the gap is large and wiow is high and far from the strong 
resonances. These considerations hold when increasing 
A'^ at a given geometry, as the localized modes do not 
change their frequency or strength of local interactions. 
The number of all resonant processes grows proportion- 
ally to the phonon density (^ N), but the strength of 
each one drops as ~ 1/ \/N for every plane- wave phonon 
involved. This leaves only the weakly coupled third-order 
resonances, near coj « wi — wiow, with a contribution 
scaling like N/^/N = y/N . In addition, there are the 



off-resonant couplings with phonons whose frequencies 
approach zero. Still, these two contributions will pose no 
problem up to at least a couple of hundreds of ions. 

Production of kink configurations in the linear trap can 
be achieved by a fast temporal variation of the transverse 
potential. Numerical simulations indicate that this pro- 
cess together with simultaneous cooling indeed leads to 
creation of kinks [Fig. [Ijb)], which remain stable at tem- 
peratures below ~ ISTdoppIci - 

Finally, our results should hold in other discrete non- 
linear models. In particular, long coherence times have 
been obtained for the FK and discrete t/)** models [sH . 

To conclude, we have demonstrated that stable clas- 
sical kink configurations exist in linear traps with > 20 
ions, as well as in circular traps. Our results suggest 
that coherence of solitonic internal modes is preserved 
for surprisingly long durations in Doppler cooled traps. 
The unique scale independent properties of the internal 
modes, specifically their gap separation from the phonon 
band, their localization to a few ions and their high fre- 
quency, suggest that such coherences can be measured 
and manipulated using existing ion trap techniques. This 
indicates that trapped-ion solitons may be useful for 
generating entanglement ll|, |32|, |33| , and implementing 



quantum information processing in large systems. 
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